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Abstract 

The asymptotic behavior of the correlator for Polyakov loop operators separated 
by a large distance R is determined for high temperature QCD. It is dominated by 
nonperturbative effects related to the exchange of magnetostatic gluons. To analyze 
the asymptotic behavior, the problem is formulated in terms of the effective field theory 
of QCD in 3 space dimensions. The Polyakov loop operator is expanded in terms 
of local gauge-invariant operators constructed out of the magnetostatic gauge field, 
with coefficients that can be calculated using resummed perturbation theory. The 
asymptotic behavior of the correlator is exp{—MR)/R, where M is the mass of the 
lowest-lying glueball in (2 -|- l)-dimensional QCD. This result implies that existing 
lattice calculations of the Polyakov loop correlator at the highest temperatures available 
do not probe the true asymptotic region in R. 



One of the basic characteristics of a plasma is the screening of electric fields. The field 
created by a static charge falls off exponentially beyond the screening radius, whose inverse is 
called the Debye mass rriD. In a quark-gluon plasma at high temperature T, chromoelectric 
fields are believed to be screened in a similar way. However it has proven to be difficult 
to give a precise definition to the Debye mass in perturbation theory. At leading order in 
the coupling constant g, rriD is proportional to gT. The next-to-leading order correction to 
mo has been calculated using a resummed perturbation theory in which gluon propagator 
corrections of order g^T"^ are summed up to all orders |l|, 0. The correction is gauge-invariant 
but infrared- divergent, indicating a sensitivity to nonperturbative effects involving the scale 

Debye screening can also be studied nonperturbatively using lattice simulations 0, 
^, . One of the simplest probes of Debye screening is the correlator of two Polyakov loop 
operators as a function of their separation R. At lowest order in resummed perturbation 
theory, the behavior of the correlator is predicted to be exp(— 2m/)i?)/i?^. By fitting the 
measured correlator to this form, one can extract a value for m^). However the resummed 
perturbation expansion is known to break down at higher orders, due to contributions that 
involve the exchange of magnetostatic gluons. These contributions, which are inherently 
nonperturbative, can be expected to dominate at distances R much greater than 1/mo- 
This raises questions about the utility of determining m/) by fitting the correlator to a 
leading order expression. 

In this paper, we determine the true asymptotic behavior of the Polyakov loop corre- 
lator at large R. We use effective field theory methods to express the Polyakov loop operator 
in terms of operators in 3-dimensional QCD. The asymptotic behavior of the Polyakov loop 
correlator is then given by a simple correlator in this effective theory. It has the form 
exp{—mHR) / R, where rriH is the mass of the lowest glueball state in (2+l)-dimensional 
QCD. The coefficient of the exponential is of order g^"^. Our result implies that the asymp- 
totic behavior of the Polyakov loop correlator is dominated by magnetostatic effects which 
have little to do with Debye screening. 

We wish to study the correlator of Polyakov loop operators in thermal QCD in 4 space- 
time dimensions with temperature T. The fundamental fields are the gauge field 74^(x, r). 
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which takes values in the SU{Nc) algebra, and the quark fields r), whose indices range 
over Nc colors and Nf fiavors. The gluon fields satisfy periodic boundary conditions in the 
Euclidean time r with period /? = 1/T, while the quark fields obey antiperiodic boundary 
conditions. The Polyakov loop operator is given by the trace of a path-ordered exponential: 




The connected part of the correlator of two Polyakov loop operators separated by a distance 
R is 

C{R) = (Lt(0) L(R)) - (L(0))l (2) 

The Polyakov loop operator (|l]) creates two or more electric gluons {Aq fields). The 
diagram for the Polyakov loop correlator (0) which is leading order in a naive expansion 
in powers of g is the 1-loop diagram in which two electrostatic gluons are exchanged. It 
gives the asymptotic behavior C{R) ~ 1/R^. This is not the correct asymptotic behavior 
for two reasons. First, thermal loop corrections generate a Debye screening mass for elec- 
trostatic gluons, so that the potential due to the exchange of an electrostatic gluon actually 
falls off exponentially, rather than like as in naive perturbation theory. Secondly, the 
true asymptotic behavior at sufficiently large temperature actually comes from higher order 
diagrams that involve the exchange of magnetostatic gluons. Our problem is to determine 
this asymptotic behavior. 

An elegant way to solve this problem is to construct a sequence of two effective field 
theories which reproduce static correlation functions at successively longer distances. Ther- 
mal QCD can be used to calculate the static correlator (Q) for any separation R. Ordinary 
perturbation theory in g"^ is accurate only for R of order 1/T or smaller, but the correlator 
can be calculated for larger R by using nonperturbative methods such as lattice gauge theory 
simulations. The first of the two effective field theories is constructed so that it reproduces 
static correlators at distances R of order 1 / {gT) or larger. In this effective theory, perturba- 
tion theory in g can be used to calculate correlators at distances of order 1 / [gT] , but lattice 
simulations are required at larger R. The second effective field theory is constructed so that it 
reproduces correlators at distances of order 1 / {g'^T) or larger. This field theory is completely 
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nonperturbative and the correlators must be calculated by lattice simulations. Nevertheless, 
we can use this field theory to determine unambiguously the asymptotic behavior of the 
Polyakov loop correlator. 

The first effective field theory, which we call electrostatic QCD (EQCD), is a 3- 
dimensional Euclidean field theory that contains the electrostatic gluon field Aq{x) and the 
magnetostatic gluon field Aj(x). Up to a normalization, they can be identified with the 
zero-frequency modes of the gluon field A^(x, r) for thermal QCD in a static gauge [0. The 
action for the effective field theory is 

5eqcd = / d^xS^hiiG.jGij) + tr(A^oA^o) + m^j tr(AoAo) + 5/:eqcd}, (3) 

where Gij = diAj — djAi + ig3[Ai,Aj] is the magnetostatic field strength, -DjAq = O^Aq + 
igs[Ai, Aq], and is the coupling constant of the 3-dimensional gauge theory. The action 
for this effective field theory is invariant under static SU{Nc) gauge transformations. If the 
fields Aq and Ai are assigned dimension 1/2, then the operators shown explicitly in (j^) have 
dimensions 3, 3, and 1. The term 5i2EQCD in (0) includes all other local gauge-invariant 
operators of dimension 2 and higher that can be constructed out of Aq and A^. The effective 
theory EQCD is completely equivalent to thermal QCD at distance scales of order l/{gT) 
or larger. The gauge coupling constant gs, the mass parameter mli, and the parameters 
in SCeqcb can be tuned as functions of T so that correlators of gauge invariant operators 
in EQCD agree with the corresponding static correlators in thermal QCD to any desired 
accuracy for R ^ 1/T By matching the two theories at tree level, we find that the 

gauge coupling constant is g^ = g{T)\/T, where g{T) is the running coupling constant at 
the momentum scale T. The mass parameter m\ in is the contribution to the square of 
the Debye screening mass from short distances of order 1/T. At leading order in g, it is 

. 2Nc + Nf n,_,_p 

ml = ^ g\T)T'. (4) 

The coefficients of some of the higher dimension operators in 5£eqcd were recently calculated 
to leading order by Chapman [§]. Since the parameters in EQCD only take into account the 
effects of the momentum scale T, they can be calculated as perturbation series in g'^{T). For 



example, the next-to-leading order correction to m^i is of order g'^{T). The Debye screening 
mass defined by the location of the pole in the gluon propagator is also given at leading 
order by (^), but mj-, has corrections of order that arise from the momentum scale gT 

The effective field theory EQCD was used by Nadkarni to study the Polyakov loop 
correlator beyond leading order in the coupling constant 0. In EQCD, the Polyakov loop 
operator is given by a simple exponential: 

L(x) = i-tr exp [-tgAo{^)/VT) . (5) 

The 1-loop diagram involving the exchange of two electrostatic gluons gives a correlator that 
falls exponentially due to electric screening: 

^^^^ ^ 8N^T- [-4^) ■ 

The corrections to this correlator of order g^ were calculated by Nadkarni [^]. They have the 
asymptotic behavior 6"^™*='-^ log(i?)/i?. In a recent reexamination of this calculation [Q, it 
has been shown that one can extract from it the correction of order g^ to the square of the 
Debye mass that was first obtained by Rebhan from the pole in the gluon propagator 

As pointed out by Nadkarni the asymptotic behavior of the Polyakov loop corre- 
lator comes not from the 1-loop diagram in which two electrostatic gluons are exchanged, 
but instead from higher-loop diagrams that involve the exchange of magnetostatic gluons. 
The simplest such diagrams are the 3-loop diagrams in Fig. 0. In perturbation theory, mag- 
netostatic gluons remain massless in EQCD, and the diagrams in Fig. |l] give a contribution 
to the correlator that falls like However this is not the correct asymptotic behavior, 

since nonperturbative effects become important at a distance R of order 1 / {g'^T) . 

In order to determine the true asymptotic behavior of the Polyakov loop correlator, it 
is useful to construct a second effective field that reproduces static correlators at distances R 
of order l/{g'^T) or larger. This effective theory, which we call magnetostatic QCD (MQCD), 
is a pure SU(3) gauge theory in 3 space dimensions. The only fields are the magnetostatic 
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gluon fields Aj(x). The action is 

Smqcb = J rf^a; l^tr(G'ijG'ij) + 5£mqcd| , (7) 

where 5£mqcd includes all local gauge-invariant operators that can be constructed out of 
Ai. The gauge coupling constant of MQCD and the parameters in 5£mqcd can be tuned so 
that MQCD is completely equivalent to EQCD, and therefore to thermal QCD, at distances 
of order l/{g^T) or larger. If g{T) is sufficiently small, the parameters of MQCD can be 
obtained by perturbative calculations in EQCD. The expansion parameter is g'^/rriei, which 
is of order g{T). The gauge theory MQCD itself is inherently nonperturbative. Any per- 
turbative expansion in powers of gs is hopelessly plagued with infrared divergences. Thus 
the correlation functions in this theory must be calculated nonperturbatively using lattice 
simulations. 

The static correlators of gauge-invariant operators in EQCD are reproduced at dis- 
tances R ^ I/ttLci by the corresponding operators in MQCD. In EQCD, the Polyakov loop 
operator creates electrostatic gluons only. It couples to magnetostatic gluons through loop 
diagrams involving electrostatic gluons, such as the 1-loop diagrams in Fig. ^. Because of 
screening, electrostatic gluons can only propagate over distances of order l/mgi. Thus for 
magnetostatic gluons with wavelengths much greater than l/rriei, the Polyakov loop opera- 
tor behaves like a point-like operator that creates magnetostatic gluons. It can therefore be 
expanded out in terms of local gauge- invariant operators constructed out of the field Ai: 

L(x) = \,{g)l + ^4^G2(x) + (8) 

where = ti (GijGij) and . . . represents operators of dimension 5 or larger, such as G^ = 
g^ti^GijGjkGki)- 

Like the parameters in the effective action (0), the coefficients in the operator expan- 
sion (H) are computable in terms of the parameters of the EQCD action using perturbation 
theory in g. Both EQCD and MQCD reproduce the nonperturbative dynamics of thermal 
QCD at distances much greater than l/rriei- Their perturbation expansions also give equiva- 
lent (although incorrect) descriptions of the long-distance dynamics. Since the coefficients in 
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the operator expansion are insensitive to the long-distance dynamics, the equivalence be- 
tween perturbative EQCD and perturbative MQCD can be exploited as a device to compute 
the coefficients. 

We proceed to calculate the coefficient Ag-z in (|^) to leading order in g. The simplest 
quantity that can be used to calculate A^a is the coupling of the Polyakov loop operator to 
two long-wavelength magnetostatic gluons, which we denote by {Q\L{0)\gg). We take the 
gluons to have momenta ki and k2, vector indices i and j, and color indices a and h. In 
perturbative MQCD, we can read off the coupling of the operator L{Q) to the two gluons 
directly from the expression (H) for the Polyakov loop operator: 

mmgg) = ^ f-ki • k25^^- + m) . (9) 

mil 

In perturbative EQCD, the coupling is given by the sum of the 1-loop diagrams in Fig. ^: 

9\ab f d^P I 1 



(0|L(0)|^^) = ^j6^'J 



(27r)3 p2 + ml (p + ki + ks)^ + ml 



(2p + ki)^(2p + 2ki + k2; 
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S'^ + ^ " r ^ ^2 2 ■ (10) 

V (p + ki)2 + m2 J 

Expanding the integrand out to second order in ki and k2 and evaluating the loop integrals, 
we find that (|10D reduces to 



Comparing (|^) and (|TT]), we can read off the coefficient Agz: 

Ag2 = (12) 
^ 3847r ^ ' 

Having determined the coefficient A^a in the operator expansion (P) to leading order 
in (yf, we can now express the Polyakov loop correlator (Q) in terms of correlators in MQCD: 

c{R) = {G\m'm) + • • • ■ (13) 



m 



The ... in ( ]I3| ) represents the contributions of higher dimension operators in the operator 
expansion (||), such as G'^(x). 

The asymptotic behavior of a correlator in MQCD, such as the one in (^3|), is related 
to the spectrum of QCD in (2+1) space-time dimensions. This is a confining gauge theory 
with a dynamically generated mass gap Mh between the vacuum and the state of next lowest 
energy. The single particle states in the spectrum are bound states of gluons (glueballs), 
and Mh is the mass of the lightest glueball. Assuming that the lowest glueball if is a scalar 
particle, the Fourier transform of the correlator (G'^(O)G'^(R)) has a pole at k"^ = —Mfj, 
as well as poles and branch cuts that are farther from the real k axis. The asymptotic 
behavior in R is dominated by the pole at k"^ = —Mfj. Denoting the residue of the pole by 
|(0|G^(0)|i7)p, the asymptotic behavior is 

f\r'2V p-M„R 

Both the mass Mh and the coupling strength {0\G'^{O)\H) can be calculated using lattice sim- 
ulations of MQCD. By dimensional analysis, Mh is proportional to g^T, while {0\G'^{O)\H) 
is proportional to [g'^T]^/'^. The overall coefficient of the exponential in ( p!^ is therefore 
proportional to g^"^. 

It is interesting to compare the asymptotic behavior (|l^ with the result one would 
obtain at leading order in perturbation theory in MQCD. The leading order diagram is 
the 1-loop diagram in which two magnetostatic gluons are exchanged. The possibility of a 
dynamically-generated magnetic screening mass mmag can be taken into account by replacing 
the propagators 1/p^ of the gluons by l/(p^ + "^mag)- The resulting expression for the 
correlator is 

C{R) ^ 2{Nl-l)\-^\ \ ^^^3 j [6 + 12m^,gi?+10(m^,gi?)2 + 4(m^,gi?)3] (15) 

The perturbative result, which is obtained by setting mmag = 0, falls like In the 

presence of magnetic screening, the asymptotic behavior is e~^™™='s^//^3 xhus this model 
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does not give the same asymptotic behavior as (p!4D, even if we make the identification 
rriH = 2mmag- 

The correlator of Polyakov loops has been calculated nonperturbatively using lattice 
simulations |, ||. Recent studies have found that, at temperatures well above the quark- 
gluon plasma phase transition, the behavior of the correlator at large R is consistent with 
the form e~^^/R^ with n in the range 1 < n < 2. At the highest temperatures available, 
the preferred value of n is close to 2, as predicted by the leading electrostatic contribution 
(^. At lower temperatures, the preferred value of n is closer to 1, consistent with the 
true asymptotic result (Q). These numerical results have a simple interpretation. The 
electrostatic contribution (^) falls off like e~^™''='^/(Ti?)^, with a coefficient proportional to 
(yf^. The magnetostatic contribution ([l^) falls off more slowly like e~^^'^^/{TR), but its 
coefficient is proportional to g^"^. At the highest temperatures that have been probed by 
lattice simulations, the running coupling constant g(T) is quite small. Because it has a 
very small coefficient, the magnetostatic contribution will probably not dominate over the 
electrostatic contribution until R is much larger than the size of present lattices. Thus the 
measured correlator behaves like e~^^/R'^, and the mass fi extracted by fitting the correlator 
to this form can be interpreted as twice the Debye screening mass. At lower temperatures, the 
magnetostatic contribution is not so strongly suppressed and the true asymptotic behavior 
Q-fj.R I ji ig probably observed on the lattice. If this interpretation of the numerical simulations 
is correct, then the mass fi extracted from fitting the correlator to the form e~^^/R has 
nothing to do with Debye screening, but instead is related to magnetostatic effects. This 
interpretation could be verified by calculating itih and (0|G^(O)|if) using a lattice simulation 
of MQCD. One could then use (JT^) to predict quantitatively how large R must be in order 
to reach the truly asymptotic region of the Polyakov loop correlator. 
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Figure Captions 

1. Three-loop diagrams in EQCD that contribute to the correlator of two Polyakov loop 
operators. The sohd lines are electrostatic gluons and the wavy lines are magnetostatic 
gluons. 

2. One-loop diagrams in EQCD that couple a Polyakov loop operator to two magneto- 
static gluons. 
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